Functional and Mechanistic Characterization of the 4,5-diepi-Isoishwarane Synthase from the Liverwort Radula lindenbergiana.

The microbial type sesquiterpene synthase RlMTPSL4 from the liverwort Radula lindenbergiana was investigated for its products, showing the formation of several sesquiterpene hydrocarbons. The main product was structurally characterized as the new compound 4,5-diepi-isoishwarane, while the side products included the known hydrocarbons germacrene A, α-selinene, eremophilene and 4,5-diepi-aristolochene. The cyclization mechanism towards 4,5-diepi-isoishwarane catalyzed by RlMTPSL4 was investigated through isotopic labeling experiments, revealing the stereochemical course for the deprotonation step to the neutral intermediate germacrene A, a reprotonation for its further cyclization, and a 1,2-hydride shift along the cascade. The absolute configuration of 4,5-diepi-isoishwarane was determined using a stereoselective deuteration approach, revealing an absolute configuration typically observed for a microbial type sesquiterpene.

Mechanistic characterisation of a sesquiterpene synthase for asterisca-1,6-diene from the liverwort Radula lindenbergiana and implications for pentalenene biosynthesis.

A sesquiterpene synthase from the liverwort Radula lindenbergiana was characterised and shown to produce the new sesquiterpene hydrocarbon (3R,9R)-asterisca-1,6-diene, besides small amounts of pentalenene. The biosynthesis of asterisca-1,6-diene was studied through isotopic labelling experiments, giving additional insights into the long discussed biosynthesis of pentalenene.

A radiation of Psylliodes flea beetles on Brassicaceae is associated with the evolution of specific detoxification enzymes.

Flea beetles of the genus Psylliodes have evolved specialized interactions with plant species belonging to several distantly related families, mainly Brassicaceae, Solanaceae, and Fagaceae. This diverse host use indicates that Psylliodes flea beetles are able to cope with different chemical defense metabolites, including glucosinolates, the characteristic defense metabolites of Brassicaceae. Here we investigated the evolution of host use and the emergence of a glucosinolate-specific detoxification mechanism in Psylliodes flea beetles. In phylogenetic analyses, Psylliodes species clustered into four major clades, three of which contained mainly species specialized on either Brassicaceae, Solanaceae, or Fagaceae. Most members of the fourth clade have broader host use, including Brassicaceae and Poaceae as major host plant families. Ancestral state reconstructions suggest that Psylliodes flea beetles were initially associated with Brassicaceae and then either shifted to Solanaceae or Fagaceae, or expanded their host repertoire to Poaceae. Despite a putative ancestral association with Brassicaceae, we found evidence that the evolution of glucosinolate-specific detoxification enzymes coincides with the radiation of Psylliodes on Brassicaceae, suggesting that these are not required for using Brassicaceae as hosts but could improve the efficiency of host use by specialized Psylliodes species.

Reinventing metabolic pathways: Independent evolution of benzoxazinoids in flowering plants.

Benzoxazinoids (BXDs) form a class of indole-derived specialized plant metabolites with broad antimicrobial and antifeedant properties. Unlike most specialized metabolites, which are typically lineage-specific, BXDs occur sporadically in a number of distantly related plant orders. This observation suggests that BXD biosynthesis arose independently numerous times in the plant kingdom. However, although decades of research in the grasses have led to the elucidation of the BXD pathway in the monocots, the biosynthesis of BXDs in eudicots is unknown. Here, we used a metabolomic and transcriptomic-guided approach, in combination with pathway reconstitution in Nicotiana benthamiana, to identify and characterize the BXD biosynthetic pathways from both Aphelandra squarrosa and Lamium galeobdolon, two phylogenetically distant eudicot species. We show that BXD biosynthesis in A. squarrosa and L. galeobdolon utilize a dual-function flavin-containing monooxygenase in place of two distinct cytochrome P450s, as is the case in the grasses. In addition, we identified evolutionarily unrelated cytochrome P450s, a 2-oxoglutarate-dependent dioxygenase, a UDP-glucosyltransferase, and a methyltransferase that were also recruited into these BXD biosynthetic pathways. Our findings constitute the discovery of BXD pathways in eudicots. Moreover, the biosynthetic enzymes of these pathways clearly demonstrate that BXDs independently arose in the plant kingdom at least three times. The heterogeneous pool of identified BXD enzymes represents a remarkable example of metabolic plasticity, in which BXDs are synthesized according to a similar chemical logic, but with an entirely different set of metabolic enzymes.

Biosynthesis, herbivore induction, and defensive role of phenylacetaldoxime glucoside.

Aldoximes are well-known metabolic precursors for plant defense compounds such as cyanogenic glycosides, glucosinolates, and volatile nitriles. They are also defenses themselves produced in response to herbivory; however, it is unclear whether aldoximes can be stored over a longer term as defense compounds and how plants protect themselves against the potential autotoxic effects of aldoximes. Here, we show that the Neotropical myrmecophyte tococa (Tococa quadrialata, recently renamed Miconia microphysca) accumulates phenylacetaldoxime glucoside (PAOx-Glc) in response to leaf herbivory. Sequence comparison, transcriptomic analysis, and heterologous expression revealed that 2 cytochrome P450 enzymes, CYP79A206 and CYP79A207, and the UDP-glucosyltransferase UGT85A123 are involved in the formation of PAOx-Glc in tococa. Another P450, CYP71E76, was shown to convert PAOx to the volatile defense compound benzyl cyanide. The formation of PAOx-Glc and PAOx in leaves is a very local response to herbivory but does not appear to be regulated by jasmonic acid signaling. In contrast to PAOx, which was only detectable during herbivory, PAOx-Glc levels remained high for at least 3 d after insect feeding. This, together with the fact that gut protein extracts of 3 insect herbivore species exhibited hydrolytic activity toward PAOx-Glc, suggests that the glucoside is a stable storage form of a defense compound that may provide rapid protection against future herbivory. Moreover, the finding that herbivory or pathogen elicitor treatment also led to the accumulation of PAOx-Glc in 3 other phylogenetically distant plant species suggests that the formation and storage of aldoxime glucosides may represent a widespread plant defense response.

Chemical profile and analysis of biosynthetic pathways and genes of volatile terpenes in Pityopsis ruthii, a rare and endangered flowering plant.

It is critical to gather biological information about rare and endangered plants to incorporate into conservation efforts. The secondary metabolism of Pityopsis ruthii, an endangered flowering plant that only occurs along limited sections of two rivers (Ocoee and Hiwassee) in Tennessee, USA was studied. Our long-term goal is to understand the mechanisms behind P. ruthii's adaptation to restricted areas in Tennessee. Here, we profiled the secondary metabolites, specifically in flowers, with a focus on terpenes, aiming to uncover the genomic and molecular basis of terpene biosynthesis in P. ruthii flowers using transcriptomic and biochemical approaches. By comparative profiling of the nonpolar portion of metabolites from various tissues, P. ruthii flowers were rich in terpenes, which included 4 monoterpenes and 10 sesquiterpenes. These terpenes were emitted from flowers as volatiles with monoterpenes and sesquiterpenes accounting for almost 68% and 32% of total emission of terpenes, respectively. These findings suggested that floral terpenes play important roles for the biology and adaptation of P. ruthii to its limited range. To investigate the biosynthesis of floral terpenes, transcriptome data for flowers were produced and analyzed. Genes involved in the terpene biosynthetic pathway were identified and their relative expressions determined. Using this approach, 67 putative terpene synthase (TPS) contigs were detected. TPSs in general are critical for terpene biosynthesis. Seven full-length TPS genes encoding putative monoterpene and sesquiterpene synthases were cloned and functionally characterized. Three catalyzed the biosynthesis of sesquiterpenes and four catalyzed the biosynthesis of monoterpenes. In conclusion, P. ruthii plants employ multiple TPS genes for the biosynthesis of a mixture of floral monoterpenes and sesquiterpenes, which probably play roles in chemical defense and attracting insect pollinators alike.

Isolation, (bio)synthetic studies and evaluation of antimicrobial properties of drimenol-type sesquiterpenes of Termitomyces fungi.

Macrotermitinae termites have farmed fungi in the genus Termitomyces as a food source for millions of years. However, the biochemical mechanisms orchestrating this mutualistic relationship are largely unknown. To deduce fungal signals and ecological patterns that relate to the stability of this symbiosis, we explored the volatile organic compound (VOC) repertoire of Termitomyces from Macrotermes natalensis colonies. Results show that mushrooms emit a VOC pattern that differs from mycelium grown in fungal gardens and laboratory cultures. The abundance of sesquiterpenoids from mushrooms allowed targeted isolation of five drimane sesquiterpenes from plate cultivations. The total synthesis of one of these, drimenol, and related drimanes assisted in structural and comparative analysis of volatile organic compounds (VOCs) and antimicrobial activity testing. Enzyme candidates putatively involved in terpene biosynthesis were heterologously expressed and while these were not involved in the biosynthesis of the complete drimane skeleton, they catalyzed the formation of two structurally related monocyclic sesquiterpenes named nectrianolins.

HDR, the last enzyme in the MEP pathway, differently regulates isoprenoid biosynthesis in two woody plants.

Dimethylallyl diphosphate (DMADP) and isopentenyl diphosphate (IDP) serves as the universal C5 precursors of isoprenoid biosynthesis in plants. These compounds are formed by the last step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway, catalyzed by (E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate reductase (HDR). In this study, we investigated the major HDR isoforms of two woody plant species, Norway spruce (Picea abies) and gray poplar (Populus × canescens), to determine how they regulate isoprenoid formation. Since each of these species has a distinct profile of isoprenoid compounds, they may require different proportions of DMADP and IDP with proportionally more IDP being needed to make larger isoprenoids. Norway spruce contained two major HDR isoforms differing in their occurrence and biochemical characteristics. PaHDR1 produced relatively more IDP than PaHDR2 and it encoding gene was expressed constitutively in leaves, likely serving to form substrate for production of carotenoids, chlorophylls, and other primary isoprenoids derived from a C20 precursor. On the other hand, Norway spruce PaHDR2 produced relatively more DMADP than PaHDR1 and its encoding gene was expressed in leaves, stems, and roots, both constitutively and after induction with the defense hormone methyl jasmonate. This second HDR enzyme likely forms a substrate for the specialized monoterpene (C10), sesquiterpene (C15), and diterpene (C20) metabolites of spruce oleoresin. Gray poplar contained only one dominant isoform (named PcHDR2) that produced relatively more DMADP and the gene of which was expressed in all organs. In leaves, where the requirement for IDP is high to make the major carotenoid and chlorophyll isoprenoids derived from C20 precursors, excess DMADP may accumulate, which could explain the high rate of isoprene (C5) emission. Our results provide new insights into the biosynthesis of isoprenoids in woody plants under conditions of differentially regulated biosynthesis of the precursors IDP and DMADP.

Predicting functions of putative fungal sesquiterpene synthase genes based on multiomics data analysis.

Sesquiterpenes (STs) are secondary metabolites, which mediate biotic interactions between different organisms. Predicting the species-specific ST repertoires can contribute to deciphering the language of communication between organisms of the same or different species. High biochemical plasticity and catalytic promiscuity of sesquiterpene synthases (STSs), however, challenge the homology-based prediction of the STS functions. Using integrated analyses of genomic, transcriptomic, volatilomic, and metabolomic data, we predict product profiles for 116 out of 146 putative STS genes identified in the genomes of 30 fungal species from different trophic groups. Our prediction method is based on the observation that STSs encoded by genes closely related phylogenetically are likely to share the initial enzymatic reactions of the ST biosynthesis pathways and, therefore, produce STs via the same reaction route. The classification by reaction routes allows to assign STs known to be emitted by a particular species to the putative STS genes from this species. Gene expression information helps to further specify these ST-to-STS assignments. Validation of the computational predictions of the STS functions using both in silico and experimental approaches shows that integrated multiomic analyses are able to correctly link cyclic STs of non-cadalane type to genes. In the process of the experimental validation, we characterized catalytic properties of several putative STS genes from the mycorrhizal fungus Laccaria bicolor. We show that the STSs encoded by the L.bicolor mycorrhiza-induced genes emit either nerolidol or α-cuprenene and α-cuparene, and discuss the possible roles of these STs in the mycorrhiza formation.

Quantifying chemodiversity considering biochemical and structural properties of compounds with the R package chemodiv.

Plants produce large numbers of phytochemical compounds affecting plant physiology and interactions with their biotic and abiotic environment. Recently, chemodiversity has attracted considerable attention as an ecologically and evolutionary meaningful way to characterize the phenotype of a mixture of phytochemical compounds. Currently used measures of phytochemical diversity, and related measures of phytochemical dissimilarity, generally do not take structural or biosynthetic properties of compounds into account. Such properties can be indicative of the compounds' function and inform about their biosynthetic (in)dependence, and should therefore be included in calculations of these measures. We introduce the R package chemodiv, which retrieves biochemical and structural properties of compounds from databases and provides functions for calculating and visualizing chemical diversity and dissimilarity for phytochemicals and other types of compounds. Our package enables calculations of diversity that takes the richness, relative abundance and - most importantly - structural and/or biosynthetic dissimilarity of compounds into account. We illustrate the use of the package with examples on simulated and real datasets. By providing the R package chemodiv for quantifying multiple aspects of chemodiversity, we hope to facilitate investigations of how chemodiversity varies across levels of biological organization, and its importance for the ecology and evolution of plants and other organisms.

Ancient origin and conserved gene function in terpene pheromone and defense evolution of stink bugs and hemipteran insects.

Insects use diverse arrays of small molecules such as metabolites of the large class of terpenes for intra- and inter-specific communication and defense. These molecules are synthesized by specialized metabolic pathways; however, the origin of enzymes involved in terpene biosynthesis and their evolution in insect genomes is still poorly understood. We addressed this question by investigating the evolution of isoprenyl diphosphate synthase (IDS)-like genes with terpene synthase (TPS) function in the family of stink bugs (Pentatomidae) within the large order of piercing-sucking Hemipteran insects. Stink bugs include species of global pest status, many of which emit structurally related 15-carbon sesquiterpenes as sex or aggregation pheromones. We provide evidence for the emergence of IDS-type TPS enzymes at the onset of pentatomid evolution over 100 million years ago, coinciding with the evolution of flowering plants. Stink bugs of different geographical origin maintain small IDS-type families with genes of conserved TPS function, which stands in contrast to the diversification of TPS genes in plants. Expanded gene mining and phylogenetic analysis in other hemipteran insects further provides evidence for an ancient emergence of IDS-like genes under presumed selection for terpene-mediated chemical interactions, and this process occurred independently from a similar evolution of IDS-type TPS genes in beetles. Our findings further suggest differences in TPS diversification in insects and plants in conjunction with different modes of gene functionalization in chemical interactions.

Biosynthesis of iridoid sex pheromones in aphids.

Iridoid monoterpenes, widely distributed in plants and insects, have many ecological functions. While the biosynthesis of iridoids has been extensively studied in plants, little is known about how insects synthesize these natural products. Here, we elucidated the biosynthesis of the iridoids cis-trans-nepetalactol and cis-trans-nepetalactone in the pea aphid Acyrthosiphon pisum (Harris), where they act as sex pheromones. The exclusive production of iridoids in hind legs of sexual female aphids allowed us to identify iridoid genes by searching for genes specifically expressed in this tissue. Biochemical characterization of candidate enzymes revealed that the iridoid pathway in aphids proceeds through the same sequence of intermediates as described for plants. The six identified aphid enzymes are unrelated to their counterparts in plants, conclusively demonstrating an independent evolution of the entire iridoid pathway in plants and insects. In contrast to the plant pathway, at least three of the aphid iridoid enzymes are likely membrane bound. We demonstrated that a lipid environment facilitates the cyclization of a reactive enol intermediate to the iridoid cyclopentanoid-pyran scaffold in vitro, suggesting that membranes are an essential component of the aphid iridoid pathway. Altogether, our discovery of this complex insect metabolic pathway establishes the genetic and biochemical basis for the formation of iridoid sex pheromones in aphids, and this discovery also serves as a foundation for understanding the convergent evolution of complex metabolic pathways between kingdoms.

Combined -omics framework reveals how ant symbionts benefit the Neotropical ant-plant Tococa quadrialata at different levels.

Ant-plant defensive mutualism is a widely studied phenomenon, where ants protect their host plants (myrmecophytes) against herbivores in return for the provision of nesting sites and food. However, few studies addressed the influence of ant colonization and herbivory on the plant's metabolism. We chose the Amazonian plant Tococa quadrialata, living in association with Azteca cf. tonduzi ants for an ant-exclusion study to reveal the chemistry behind this symbiosis. We found that colonized plants did not only benefit from protection but also from increased amino acid and nitrogen content, enabling better performance even in an herbivore-free environment. In contrast, ant-deprived T. quadrialata plants accumulated more ellagitannins, a major class of constitutive defense compounds. Moreover, herbivory-induced jasmonate-mediated defense responses, including the upregulation of signaling and defense genes and the emission of volatiles irrespective of colonization status. Altogether, we show how ant-colonization can influence the general and defense-related metabolism and performance of myrmecophytes.

Augmented Enterocyte Damage During Candida albicans and Proteus mirabilis Coinfection.

The human gut acts as the main reservoir of microbes and a relevant source of life-threatening infections, especially in immunocompromised patients. There, the opportunistic fungal pathogen Candida albicans adapts to the host environment and additionally interacts with residing bacteria. We investigated fungal-bacterial interactions by coinfecting enterocytes with the yeast Candida albicans and the Gram-negative bacterium Proteus mirabilis resulting in enhanced host cell damage. This synergistic effect was conserved across different P. mirabilis isolates and occurred also with non-albicans Candida species and C. albicans mutants defective in filamentation or candidalysin production. Using bacterial deletion mutants, we identified the P. mirabilis hemolysin HpmA to be the key effector for host cell destruction. Spatially separated coinfections demonstrated that synergism between Candida and Proteus is induced by contact, but also by soluble factors. Specifically, we identified Candida-mediated glucose consumption and farnesol production as potential triggers for Proteus virulence. In summary, our study demonstrates that coinfection of enterocytes with C. albicans and P. mirabilis can result in increased host cell damage which is mediated by bacterial virulence factors as a result of fungal niche modification via nutrient consumption and production of soluble factors. This supports the notion that certain fungal-bacterial combinations have the potential to result in enhanced virulence in niches such as the gut and might therefore promote translocation and dissemination.

Dynamic regulation of volatile terpenoid production and emission from Chrysanthemum morifolium capitula.

Flower-associated communities consist of both mutualistic and antagonistic organisms. We have limited knowledge on how flowers regulate volatiles to balance their defense against antagonists and the attraction of beneficial organisms necessary for reproductive success. Asteraceae is the largest family among flowering plants. Its representatives are characterized by unique inflorescence called capitulum, which has been reduced to a reproduction unit resembling a single flower. Here, we chose Chrysanthemum morifolium, a model species of Asteraceae, to investigate how the capitulum balances the accumulation and emission of floral terpenoid volatiles that are implicated in defense and pollinator attraction, respectively. Our results showed that the capitula of C. morifolium produce and emit complex mixtures of monoterpenoids and sesquiterpenoids. The highest concentrations of terpenoids were detected in the bud stage of the capitula. In contrast, the capitulum reached the highest emission level prior to full blooming. The disc florets were the dominant organs of terpenoid accumulation and emission in the full-openness stage. To understand the molecular basis of volatile terpenoid biosynthesis in C. morifolium, experiments were designed to study terpene synthase (TPS) genes, which are pivotal for terpene biosynthesis. Eight CmCJTPS genes were identified in the transcriptomes of C. morifolium, and the proteins encoded by five genes were found to be biochemically functional. CmCJTPS5 and CmCJTPS8 were the multi-product enzymes catalyzing the monoterpenoid and sesquiterpenoid formation, which closely matched the major terpenoids produced in the flower heads. The five functional terpene synthase genes exhibited similar temporal expression patterns but diverse spatial expression levels, suggesting tissue-specific functions. Altogether, our results illustrate the dynamic patterns of accumulation and emission of floral volatile terpenoids implicated in defense and attracting pollinators in C. morifolium, for which both the regulation of TPS gene expression and the regulation of release may play critical roles.

Origin and early evolution of the plant terpene synthase family.

As a midsized gene family conserved more by lineage than function, the typical plant terpene synthases (TPSs) could be a valuable tool to examine plant evolution. TPSs are pivotal in biosynthesis of gibberellins and related phytohormones as well as in formation of the extensive arsenal of specialized plant metabolites mediating ecological interactions whose production is often lineage specific. Yet the origin and early evolution of the TPS family is not well understood. Systematic analysis of an array of transcriptomes and sequenced genomes indicated that the TPS family originated after the divergence of land plants from charophytic algae. Phylogenetic and biochemical analyses support the hypothesis that the ancestral TPS gene encoded a bifunctional class I and II diterpene synthase producing the ent-kaurene required for phytohormone production in all extant lineages of land plants. Moreover, the ancestral TPS gene likely underwent duplication at least twice early in land plant evolution. Together these two gave rise to three TPS lineages leading to the extant TPS-c, TPS-e/f, and the remaining TPS (h/d/a/b/g) subfamilies, with the latter dedicated to secondary rather than primary metabolism while the former two contain those genes involved in ent-kaurene production. Nevertheless, parallel evolution from the ent-kaurene­producing class I and class II diterpene synthases has led to roles for TPS-e/f and -c subfamily members in secondary metabolism as well. These results clarify TPS evolutionary history and provide context for the role of these genes in producing the vast diversity of terpenoid natural products observed today in various land plant lineages.

The wheat dioxygenase BX6 is involved in the formation of benzoxazinoids in planta and contributes to plant defense against insect herbivores.

Benzoxazinoids are plant specialized metabolites with defense properties, highly abundant in wheat (Triticum), one of the world's most important crops. The goal of our study was to characterize dioxygenase BX6 genes in tetraploid and hexaploid wheat genotypes and to elucidate their effects on defense against herbivores. Phylogenetic analysis revealed four BX6 genes in the hexaploid wheat T. aestivum, but only one ortholog was found in the tetraploid (T. turgidum) wild emmer wheat and the cultivated durum wheat. Transcriptome sequencing of durum wheat plants, damaged by either aphids or caterpillars, revealed that several BX genes, including TtBX6, were upregulated upon caterpillar feeding, relative to the undamaged control plants. A virus-induced gene silencing approach was used to reduce the expression of BX6 in T. aestivum plants, which exhibited both reduced transcript levels and reduced accumulation of different benzoxazinoids. To elucidate the effect of BX6 on plant defense, bioassays with different herbivores feeding on BX6-silenced leaves were conducted. The results showed that plants with silenced BX6 were more susceptible to aphids and the two-spotted spider mite than the control. Overall, our study indicates that wheat BX6 is involved in benzoxazinoid formation in planta and contributes to plant resistance against insect herbivores.

Evolution of DIMBOA-Glc O-Methyltransferases from Flavonoid O-Methyltransferases in the Grasses.

O-Methylated benzoxazinoids (BXs) and flavonoids are widespread defenses against herbivores and pathogens in the grasses (Poaceae). Recently, two flavonoid O-methyltransferases (FOMTs), ZmFOMT2 and ZmFOMT3, have been reported to produce phytoalexins in maize (Zea mays). ZmFOMT2 and ZmFOMT3 are closely related to the BX O-methyltransferases (OMTs) ZmBX10-12 and ZmBX14, suggesting a common evolutionary origin in the Poaceae. Here, we studied the evolution and enzymatic requirements of flavonoid and BX O-methylation activities in more detail. Using BLAST searches and phylogenetic analyses, we identified enzymes homologous to ZmFOMT2 and ZmFOMT3, ZmBX10-12, and ZmBX14 in several grasses, with the most closely related candidates found almost exclusively in species of the Panicoideae subfamily. Biochemical characterization of candidate enzymes from sorghum (Sorghum bicolor), sugar cane (Saccharum spp.), and teosinte (Zea nicaraguensis) revealed either flavonoid 5-O-methylation activity or DIMBOA-Glc 4-O-methylation activity. However, DIMBOA-Glc 4-OMTs from maize and teosinte also accepted flavonols as substrates and converted them to 3-O-methylated derivatives, suggesting an evolutionary relationship between these two activities. Homology modeling, sequence comparisons, and site-directed mutagenesis led to the identification of active site residues crucial for FOMT and BX OMT activity. However, the full conversion of ZmFOMT2 activity into BX OMT activity by switching these residues was not successful. Only trace O-methylation of BXs was observed, indicating that amino acids outside the active site cavity are also involved in determining the different substrate specificities. Altogether, the results of our study suggest that BX OMTs have evolved from the ubiquitous FOMTs in the PACMAD clade of the grasses through a complex series of amino acid changes.

Comparative Genomic and Metabolomic Analysis of Termitomyces Species Provides Insights into the Terpenome of the Fungal Cultivar and the Characteristic Odor of the Fungus Garden of Macrotermes natalensis Termites.

Macrotermitinae termites have domesticated fungi of the genus Termitomyces as food for their colony, analogously to human farmers growing crops. Termites propagate the fungus by continuously blending foraged and predigested plant material with fungal mycelium and spores (fungus comb) within designated subterranean chambers. To test the hypothesis that the obligate fungal symbiont emits specific volatiles (odor) to orchestrate its life cycle and symbiotic relations, we determined the typical volatile emission of fungus comb biomass and Termitomyces nodules, revealing α-pinene, camphene, and d-limonene as the most abundant terpenes. Genome mining of Termitomyces followed by gene expression studies and phylogenetic analysis of putative enzymes related to secondary metabolite production encoded by the genomes uncovered a conserved and specific biosynthetic repertoire across strains. Finally, we proved by heterologous expression and in vitro enzymatic assays that a highly expressed gene sequence encodes a rare bifunctional mono-/sesquiterpene cyclase able to produce the abundant comb volatiles camphene and d-limonene. IMPORTANCE The symbiosis between macrotermitinae termites and Termitomyces is obligate for both partners and is one of the most important contributors to biomass conversion in the Old World tropic's ecosystems. To date, research efforts have dominantly focused on acquiring a better understanding of the degradative capabilities of Termitomyces to sustain the obligate nutritional symbiosis, but our knowledge of the small-molecule repertoire of the fungal cultivar mediating interspecies and interkingdom interactions has remained fragmented. Our omics-driven chemical, genomic, and phylogenetic study provides new insights into the volatilome and biosynthetic capabilities of the evolutionarily conserved fungal genus Termitomyces, which allows matching metabolites to genes and enzymes and, thus, opens a new source of unique and rare enzymatic transformations.